Effective interactions in colloid - semipermeable membrane systems 



in 
o 
o 

(N 
o 

Q 

(N 



o 



c 

o 
o 



> 

On 
(N 

in 
o 

-I— > 



i 

o 
o 



x 

S3 



Pawel Bryk 

Department for the Modeling of Physico-Chemical Processes, 
Maria Curie- Sklodowska University, 20-031 Lublin, Polan^j 
(Dated: February 2, 2008) 

We investigate effective interactions between a colloidal particle, immersed in a binary mixture 
of smaller spheres, and a semipermeable membrane. The colloid is modeled as a big hard sphere 
and the membrane is represented as an infinitely thin surface which is fully permeable to one of the 
smaller spheres and impermeable to the other one. Within the framework of the density functional 
theory we evaluate the depletion potentials, and we consider two different approximate theories - 
the simple Asakura-Oosawa approximation and the accurate White-Bear version of the fundamen- 
tal measure theory. The effective potentials are compared with the corresponding potentials for a 
hard, nonpermeable wall. Using statistical-mechanical sum rules we argue that the contact value 
of the depletion potential between a colloid and a semipermeable membrane is smaller in magni- 
tude than the potential between a colloid and a hard wall. Explicit calculations confirm that the 
colloid-semipermeable membrane effective interactions are generally weaker than these near a hard 
nonpermeable wall. This effect is more pronounced for smaller osmotic pressures. The depletion 
potential for a colloidal particle inside a semipermeable vesicle is stronger than the potential for the 
colloidal particle located outside of a vesicle. The asymptotic decay of the depletion potential for the 
semipermeable membrane is similar to that for the nonpermeable wall and reflects the asymptotics 
of the total correlation function of the corresponding binary mixture of smaller spheres. Our results 
demonstrate that the ability of the membrane to change its shape constitutes an important factor 
in determining the effective interactions between the semipermeable membrane and the colloidal 
macroparticle. 



I. INTRODUCTION 

Depletion interactions play an important role in deter- 
mining the stability of colloidal systems^. They arise 
in asymmetric mixtures 3 due to a tendency to maximize 
volume available to centers of smaller species. In suit- 
ably prepared colloidal suspensions, such that only bare 
hard-core interactions are left, these effective interactions 
are attractive at small distances and can exhibit a poten- 
tial barrier and an oscillatory tail for large separations^. 
Depletion forces also arise if a big colloidal particle sus- 
pended in a sea of small particles approaches a fixed ob- 
ject such as a planar wall. 

The first theoretical approach to calculate the effective 
interactions was performed by Asakura and Oosawa£, 
and later, by Vrijt The resulting effective potential 
is attractive at small distances and proportional to the 
packing fraction of the depleting agent. More recent 
approaches include various integral equation&Siifiiii and 
density functional theoriesi2*i2ii4. Several experimen- 
tal techniques have been developed to investigate deple- 
tion interactions. Total internal reflection microscopy 
has been used to investigate depletion potentials be- 
tween a wall and a sphere immersed in a sea of spherical 
polymersi^ and rod-like particleisi&— Video-microscopy 
has been used to measure the depletion potential between 
two big particles^, the potential of a big colloid near 
an edge of a terraced and the potential of a big colloid 
sphere inside a vesicle— In the latter experiment it has 
been observed that the big particle is attracted to regions 
of high local curvature. This attraction may in turn in- 
duce transformation of the shape of a membrane. Thus 
the depletion interactions between a colloidal particle and 



a membrane may be helpful in a better understanding of 
many biological processes. 

Theoretical attempts to describe the colloid- vesicle ex- 
periment were initiated by Roth et alM who used den- 
sity functional theory (DFT) to study depletion forces 
between a big hard sphere, immersed in a sea of smaller 
spheres, and a hard, nonpermeable spherical wall. These 
authors also presented a simple framework for tackling 
cavities of a general shape. A more refined approach 
was presented very recently 23,24 . Bickel investigated de- 
pletion interactions near soft, fluctuating surfaces and 
showed a mechanism for an encapsulation of the colloidal 
macroparticle 2 ^. In all these studies the membrane itself 
is assumed to be nonpermeable. However, an inherent 
attribute of a semipermeable membrane system is the 
presence of a fluid on both membrane sides. The dif- 
ference in composition of the fluid on both sides, which 
appears due to the selective permeability of a membrane, 
gives rise to the osmotic pressure. 

Several microscopic models of semipermeable mem- 
brane systems have been proposed in the literature— 
The simplest model is the so-called ideal membrane - an 
infinitely thin surface which is fully permeable to some of 
the species and nonpermeable to others. The fluid struc- 
ture near such a membrane can be ascertained by means 
of firs1i2L2& and second order— integral equation theo- 
ries or DFT—. A somewhat more sophisticated model 
accounts for finite thickness by treating the membrane 
as a potential barrier of a finite height and width 31-32 33 . 
Finally, in some studies the discrete structure of a mem- 
brane was also taken into account by considering mem- 
branes built of particles arranged either randomly 3 ^ or in 
a regular manner— ii&SiiSI. 
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In the present paper we investigate the effective in- 
teractions between a semipermeable membrane and a 
big colloidal macroparticle. We use the reliable DFT 
approach that was used previously to study depletion 
potentials near impenetrable surfaces and inside hard- 
sphere cavities. This enables us to make a direct com- 
parison and elucidate the effect of the semipermeability 
on the effective interactions. The paper is arranged as 
follows: in Section^ we consider sum rules for semiper- 
meable membranes. Section UTTl outlines the theories ap- 
plied in the present paper. Section llVI presents results 
and we conclude in Section Ivl 



II. SUM RULES FOR SEMIPERMEABLE 
MEMBRANES 




Consider two component mixture of species A and B 
and diameters a a and as, respectively. The membrane is 
fully permeable to the component A but is nonpermeable 
to the component B. In the course of the derivation we 
consider the ideal membrane of the spherical shape - a 
simple model for a semipermeable vesicle. The sum rule 
for semipermeable planar membrane is then obtained by 
taking the limit Rm — > oo. The membrane is modeled 
by imposing an external potential of the form 

Vr*(r) S 0; VTCr) = { ^ ^Al . 0) 

In Eq. JU Rm — Rm + ct -b/2, and denotes the ac- 
tual membrane radius (the dashed circle in Fig. QJ. Rm 
stands for the radius at which the profile of the nonper- 
meable species makes a jump (the solid circle in Fig.[Q. 
The membrane separates the two parts of the system: 
part II with a two-component fluid and part I, available 
only to the permeable component. The bulk densities in 
part II, p^ 1 ^ and Pg 1 ^, are the input to the calculations. 
The system is in the termodynamical but not mechanical 
equilibrium. The pressure p is different in the two parts 
of the system and we define the osmotic pressure as 

n = P (»Hti I \ P ^)- P ^). (2) 

The bulk density in the one-component part of the sys- 
tem, Pa, is obtained by imposing the equality of the 
chemical potentials of the permeable component in both 
parts of the system 

P { l I \pT\pf ) ) = ti\ttY (3) 

It is clear that the osmotic pressure defined by Eqs. 
and J2J depends only on the equation of state of the 
system and is independent of the shape of the membrane. 

We decompose the grand canonical potential of the sys- 
tem SI into a reference, bulk part Q bulk = — pC-OyC- 17 ) — 
pCOyW, where V^ 1 ' and V^- 11 ' are the volumes of parts 
I and II, respectively, and to an excess part Vt ex : SI = 



FIG. 1: Schematic plot of a semipermeable vesicle system. 

Qbuik _|_ Qex _ \y e choose Rm and the corresponding 
spherical surface of the area &kR m as the dividing in- 
terface. This implies that the volumes of both parts are: 
= 4/37ri?|f, and = 4/3ir{7t - R 3 M ), whereas 

the densities in the reference system are 

Pa lk = P ( I ] Q(Rm -r) + p A a) e(r - R M ) 
p b B ulk = p l ^Q(r-R M ), (4) 

where 9 is the Heaviside unit-step function (cf. dotted 
lines in Fig.^l. R is a macroscopically large radius con- 
sidered in the thermodynamical limit. Denoting the ex- 
cess grand potential per unit are as the surface tension 
7, we have 

n = 1 {r m )^R 2 m -p {u) v^ u \Rm)-p^v^{Rm) . 

(5) 

We derive the sum rules for semipermeable membranes 
by expressing the grand potential as a functional of the 
local densities 




(6) 

where F is the Hclmholtz free energy functional that will 
be defined later. The only explicit dependence of SI on 
Rm is contained in the second term in the right-hand side 
of Eq. ©. At thermodynamical equilibrium the grand 
potential satisfies the Euler-Lagrange equation 

Sn[p A ,p B ] = 5Q.[pa,Pb] = „ (7] 
S PA (r) 5 PB (r) ■ 1 ' 

The sum rules are obtained by taking the derivative 
of the grand potential Eq. @) with respect to Rm, at 
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constant chemical potentials fii and temperature T 

( dn \ r 3 gOl^pg] dpi(r) 

\9Rm)^ t JfJ T S Pl (r) dR M 

: A d ** 



i=A,B ' 



(8) 



By virtue of Eq. J7J) the first sum on the right-hand 
side of Eq. (JSJ vanishes whereas the second term yields 
0~ 1 4:'kR\ i Pb (R m ), where pb{R m ) * s the contact density 
of the nonpermeable component of a fluid mixture at the 
ideal semipermeable vesicle. Thus the sum rule reads 



on 

8RmJ 



) = 4nR M p(R M 



), 



(9) 



and is satisfied by all density functionals within weighted- 
density approximation22iiMi. 

Taking the derivative of Eq. with respect to Rm 
and inserting to Eq. @ the sum rule can be expressed as 



Pb(R m ) 



fm. 



VHRm) 
Rm 



+ (3 
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(10) 



In the limit Rm — * oo only the first term of Eq. i|10|) 
survives and we recover the planar resuhS&iS^Ma that 
the contact value of the density profile of the nonperme- 
able species at the planar semipermeable membrane is 
equal to the osmotic pressure. Since II is always lower 
than p 1 - 11 ^ we conclude that for the same fluid the con- 
tact value at the membrane is lower than the sum of the 
contact values at a nonpermeable hard wall 45 . One can 
also consider the case in which the two component mix- 
ture is inside the vesicle and the permeable component 
is present outside the vesicle, i.e. the "reverse" of the 
situation plotted in Fig. ^ In this case the contact value 
of the nonpermeable component will be higher than that 
at the planar membrane. 



III. THEORY FOR DEPLETION POTENTIALS 

We consider a single big colloidal hard-sphere of diam- 
eter ac immersed in a two-component mixture of smaller 
hard-sphere of diameters a a and ob, respectively. The 
system is in contact with an ideal semipermeable mem- 
brane which is permeable to the A component but non- 
permeable to the B component, as well as to the big 
sphere C. We also assume that a a < <?b < &c- The 
goal is to evaluate the effective depletion potential W 
between the big colloidal hard-sphere and the semiper- 
meable membrane. In our study we use a versatile ap- 
proach proposed in Refsii^fi^. According to this theory 
the depletion potential is evaluated on the basis of the 
density profiles of the small spheres undisturbed by the 
presence of the big sphere. In the dilute limit of the big 
spheres (pc — ► or, equivalently, \xc — ► — oo) the deple- 
tion potential IF at a point r can be expressed in terms of 



the difference between the one-particle direct correlation 
function Cq in the bulk (r — » oo) and at r 

/3W(r)= lim ( cL X) (r -f oo) - <j$ (r)) . (11) 

A convenient route for the evaluation of Cq (r) is pro- 
vided by density functional theory^. The focal point 
in this approach is the excess (over the ideal gas) free 
energy, F ex = F - £ J d 3 r{p l (r)(ln( ( o 4 (r)) - 1)}. The 
one-particle direct correlation function is then directly 
accessible via 



cg>(r) = -/J 



5F e: 



5p c {r) 



(12) 



For mixtures of hard-spheres the exact excess free en- 
ergy formula is not known and some approximations are 
required. In this work we consider two approximate func- 
tionals. The first one is the low density limit functional 

P F «° = 4E [ d3r I dh'p^p^f^v-v') , (13) 

where /y is the Mayer bond between a particle of species 
i and a particle of species j. In the low density limit the 
density profiles Pa( t ) and Pb(?) tend to p^ and p h g lk , 
respectively. Then from Eqs. I|ll(l - I|13|l we obtain the 
Asakura-Oosawa approximation for the depletion poten- 
tial 



(3W AU (r) 



£ / dV(p^- Pj (//) )/,c(r-r') 

=A,B J 



(14) 

where /ic(r — r') = — 8((a i + ac)/2 — |r — r'|). Equation 
(|14fl implies that W AO is identically zero for distances 
larger than <jb from the membrane. 

The second approximate functional considered in this 
work is Rosenfeld's Fundamental Measure Theory^ 
(FMT). Within this approach the excess free energy of 
an N component hard-sphere mixture is given by 



l3F ex = / d 3 r $({n Q }) , 



where n a denote weighted densities 

N 
3=1 ' 



N 

*«(*)= E J d^'pA*') w${i 



(15) 



(16) 



The weight functions Wa depend of geometrical proper- 
ties of individual species, see earlier papers for explicit 
formulas^Z^. 

From a number of prescriptions for the excess free en- 
ergy density for the given problem we have chosen the 
elegant and inspiring White-Bear (WB) formula™^ 



H{n a }) = -ra ln(l-n 3 ) + 



+{n\ - 3n 2 n V2 ■ n V2 ) 



nin 2 - nyi • n V2 
1 — rz-3 

n 3 + (l-n 3 ) 2 ln(l-n 3 ) 
367r(n 3 ) 2 (l-n 3 ) 2 



(17) 



4 



0.05 



. 1 1 1 1 1 


' 1 ' 1 ' 1 ' - 
(a)_ 

- 

\ A 


, I.I, 


v ■ ■ 


, 1.1,1. 



10 



12 



14 



16 18 20 22 

r/a. 



24 



0.11 



V o.i 



0.09 - 




FIG. 2: (a) Density profiles of a hard-sphere mixture near a 
semipermeable vesicle. The bulk densities are p^^a A = 0.1 
and pg 1 oa = 0.05. (b) The contact value of the density 
profile of the nonpermeable component (solid line) and the 
right-hand side of the sum rule, Ea llOl (circles) plotted as a 
function of the inverse membrane radius, Rm. 



As explained ini&ll FMT provides a very nice method 
of calculating the depletion potentials. In the present 
case of a ternary hard-sphere mixture the dilute limit 
of species C can be obtained by considering the set of 
restricted weighted densities 



tin 



(r) = £ fd 3 r' Pj (r')w^(r-v') 

„■ A D ** 



(18) 



j=A,B' 

The depletion potential is then obtained from 



(3W WB (r) 



where 



Ht a (i>; {<})=/? 



d\^ a {v'-{nf })wP(r-r') , (19) 



0$({n°"}) 
dnt il 



^(tf}) 
dnt il 



(20) 

The computational strategy consists of two stages. First 
the profiles for the binary A+B mixture are evaluated 
by minimizing the White-Bear version of the fundamen- 
tal measure functional. The resulting set of restricted 
weighted densities is then used as an input to a sepa- 
rate program, in which the depletion potentials are com- 
puted by applying Eqs. {T5J-IJ2DJ. The strategy out- 
lined above offers an important advantage, because the 
computational problem is effectively reduced from three- 
component to two-component mixture. This method is 
especially appreciated when the density profile is a func- 
tion of more than one scalar—. 
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FIG. 3: Depletion potentials for a single big hard sphere 
C immersed in a mixture of hard spheres A+B of diameters 
o~b = 2a a- The diameter of the big hard sphere ac = 3cta- 
The solid lines denote depletion potentials between a big hard 
sphere and a semipermeable planar membrane evaluated from 
the WB theory. The dashed lines denote depletion potentials 
for the same mixture but at a hard nonpenetrable wall evalu- 
ated from the WB theory. The dotted lines denote depletion 
potentials at a semipermeable planar membrane evaluated 
using AO approximation. Packing fractions are r\A = 0.05, 
r\B = 0.35 in (a) and t]a = 0.15, 77s = 0.25 in (b). 



IV. RESULTS AND DISCUSSION 

A. Depletion Potentials 

Typical density profiles resulting from the WB the- 
ory are displayed in the upper panel of Fig. [3 The 
profiles were evaluated for a binary hard sphere mix- 
ture of sizes o~b = 2a a near a semipermeable vesi- 
cle of the size Rm = I60V1. The bulk densities are 

0.05. The osmotic pressure 



(II) 3 
Pa °a 



0.1 and p% a\ 



PUa A =0.105 391 59. The contact value of the density 
profile pb{R\[) =0.102 373 22 whereas the right hand- 
side of Eq. (Unj is 0.102 373 78. Similar level of agreement 
was found for other radii and this confirms high accuracy 
of the numerical results. The plot of the contact value as 
a function of the inverse membrane radius evaluated for 
these bulk densities is shown in the lower panel of Fig.|3 
In the planar membrane limit, aA/Ru = the contact 
value tends to the osmotic pressure. 

We turn now to the depletion potentials. In order to 
restrict a rather large parameter space we consider the 
case of the fixed diameters of smaller spheres = la a- 
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FIG. 4: The same as in Fig.|3]but now for packing fractions 
t\a = 0.25, rjB = 0.15 in (a) and r\A = 0.35, r\B = 0.05 in (b). 
The insets show the zoom-out of the main figures. 
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FIG. 5: Depletion potentials for a single big colloidal sphere 
of diameter ac = 5<ja near a semipermeable planar mem- 
brane. The packing fraction of the nonpermeable component 
r\B = 0.15, and o\b = la A- The solid lines denote the po- 
tentials for the packing fractions of the permeable component 
TjA = 0.10, 0.15, 0.20, and 0.25. The arrow indicates the in- 
crease of r\A- The dashed line indicates the result for t\a — 0. 



Moreover, we consider the case of a constant total pack- 
ing fraction rj = Y,i=A B ^ = Ei=A,s (V^Pi '°f = °- 4 
and the case of a fixed t\b = 0.15. By varying packing 
fraction of the components we have a possibility to in- 
vestigate fairly dense systems but with different osmotic 
pressures. Thus we hope to draw rather general conclu- 
sions. 

Figures El and ^ show the depletion potentials between 
a big colloidal particle of the size ac — 3(7 a and a planar 
semipermeable membrane, plotted as a function of the 
distance to the membrane surface z. The potentials were 
evaluated for the constant total packing fraction rj = 0.4. 

The solid lines are the results of the WB theory 
whereas the dotted lines denote the results of the AO ap- 
proximation. The dashed lines denote W WB but for the 
planar, nonpenetrable wall. We observe that the differ- 
ences between the colloid-wall and the colloid-membrane 
potentials are rather significant. In all cases the effect 
of the membrane permeability is to decrease the effective 
interactions, relative to the hard wall. The oscillations 
are less pronounced and the contact values of the deple- 
tions potentials are smaller in magnitude. These trends 
are already well visible for small packing fractions of the 
smaller spheres (cf. Fig. f]A = 0.05, rjs = 0.35) and 
increase with tja (cf. Fig. 03, t}a = 0.15, tjb = 0.25). For 
high packing fractions of the smaller spheres and low tjb, 
(see Fig. QJ the difference between the corresponding po- 
tentials is rather dramatic. The contact value of (3W for 
the membrane is less than -0.5 whereas for the hard wall 
this value is around -4 (see the insets to Fig. 0J . Note 



also that the oscillations, and in particular the first max- 
imum, are significantly smaller for the membrane than 
for the wall. 

That the depletion potential near a semipermeable 
membrane is weaker than the potential near a nonper- 
meable wall can be anticipated already from Eq. I|10(l . 
For the one-component hard sphere system the contact 
value of the depletion potential, as predicted by the ac- 
curate parameterization^, is linear in the packing frac- 
tion of the small spheres. The packing fraction of the 
small spheres itself is related to the contact value of the 
profile near the wall, thus the higher the pressure and 
the contact value of the density profile, the stronger the 
depletion potential. While for mixtures such parameter- 
ization is unknown it is reasonable to assume a similar 
behavior. Since the contact value near the membrane is 
always lower, the depletion potential is expected to be 
weaker and this is confirmed by the explicit calculations. 

It is interesting to note that the simple Asakura- 
Oosawa approximation (dotted lines if Figs. and 0J 
performs relatively well when compared with the WB 
theory. It is clear that the AO approximation cannot ac- 
count for the oscillatory part of the depletion potential 
but the contact value resulting from this approach does 
not differ too much from the accurate WB theory. 

Figure [5] shows the depletion potentials between a sin- 
gle big colloidal particle of the size ac — 5aA and a 
semipermeable planar membrane, evaluated for a con- 
stant packing fraction of the nonpermeable component 



6 




FIG. 6: (a) Depletion potentials between a single big colloidal 
sphere of diameter ac = 3<ta and a semipermeable vesicle of 
the radius Rm = 20cta plotted as a function of the relative 
distance to the vesicle (r — Rm)/ga- The solid line denotes 
the potential for the system in which the two component part 
of the system is located outside of the vesicle whereas the 
dashed line denotes the potential the system in which the two 
component part of the system is located inside the vesicle. 
The packing fractions r\A = 0.05 and t\b = 0.35. (b) Contact 
value of the depletion potential between a single big colloidal 
sphere of diameter ac = 5<ta and a semipermeable vesicle 
plotted as a function of the vesicle curvature. The solid line 
denotes the potential for the system in which the two com- 
ponent part of the system is located outside of the vesicle 
whereas the dashed line denotes the potential the system in 
which the two component part of the system is located inside 
the vesicle. The packing fractions t\a = 0.30 and t\b = 0.10. 
Black dot indicates the planar result limit. 



r/B = 0.15 and for different packing fractions of the per- 
meable component t]a — 0.10, 0.15, 0.20, and 0.25. The 
direction of the increase of t\a is indicated by the arrow 
and it is also the direction of decreasing osmotic pressure. 
We observe that for constant r\B the effective interactions 
get stronger with an osmotic pressure increase (the same 
trend is observed in Figs. l3l4|) . For reference we also show 
the potential for tja — (dashed line) and it is clearly 
the strongest potential for this tjb- 

Finally in Fig. HJ we show the influence of the curva- 
ture of the membrane on the depletion potentials. We 
consider the depletion potential near a convex and a con- 
cave semipermeable membrane. This corresponds to the 
cases where the two component part of the system is 
located outside of, and inside the vesicle, respectively. 
Similar to the hard wall case22 we find that the deple- 
tion potential inside the vesicle is stronger, whereas the 



potential outside of the vesicle weaker, than the corre- 
sponding potential at the planar membrane (cf. Fig- Et)- 
However, this effect is moderated in comparison to the 
hard wall case because, in general, the depletion poten- 
tials near semipermeable membranes are weaker. Figure 
Fig. 03 shows the contact value of the depletion poten- 
tial at a semipermeable membrane plotted as a function 
of the inverse vesicle radius. We find that, although the 
convergence to the planar membrane result (marked by 
the black dot) is rather slow, one has to go to rather small 
radii in order to noticeably increase/decrease the contact 
value of the depletion potential. 



B. Asymptotic behavior 

The behavior of the depletion potential at large sep- 
arations can be ascertained by using the general theory 
of asymptotics of pair correlations functions^S. It was 
demonstrated that the decay of many structural proper- 
ties of inhomogeneous fluids such as density profile, sol- 
vation force and depletion potential^ should be the same 
as the asymptotics of the bulk pair correlation functions 



Recently Grodon et 



i 54,55 



presented a detailed in- 



vestigation of the asymptotic decay of pair correlation 
functions in binary hard-sphere mixtures. The main find- 
ing is that for sufficiently asymmetric mixtures there is 
a structural crossover line in the (tiaMb) plane which 
separates two distinct regimes of the asymptotic decay. 
Below we recall main predictions of the general theory of 
asymptotics of pair correlations functions for mixtures. 

The Fourier transforms of the total correlation func- 
tions for a binary mixture of A and B species can be 
written as 



hij(k) 



D(k) 



i,j = A,B . 



(21) 



While the nominators Nij(k) in Eq. (|21|l are different 
for different pairs ij, all three total correlation functions 
have the same denominator D(k) 



PB&B B {k)] 



PAPBC A ' B {k) . 

(22) 

In the above c\j (k) is the Fourier transform of the (two- 
particle) direct correlation function. The inverse Fourier 
transform of Eq. (|22|l is formally obtained via the residue 
theorem. While in general an infinite number of residues 
and the corresponding poles contributes to hij(r), the 
asymptotic decay is governed by the leading order pole 
p(LO) _ _j_ ai _|_ j aQ ^ that is by the pole with the smallest 
imaginary part. Close to a crossover point two poles 
j/ 1 ) = a\ +iao and p' 2 ' = a\ + iclq have similar imaginary 
parts and both contribute to the decay. Provided that the 
Fourier transforms of cy (r) are known analytically, the 
poles can be conveniently determined by finding complex 
solutions to the equation D(k) = 0. In the present case 
of the binary hard-sphere mixture the direct correlation 
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FIG. 7: Asymptotic decay of the depletion potential between 
a single big colloidal particle and a semipermeable planar 
membrane. Shown are the results of the WB theory (solid 
lines) and the asymptotic forms evaluated from Eqs. 1231 - 
(12 It (dashed lines), oc = 3<r a and the packing fractions are 
r\A = 0.05 and r\s = 0.35 for the potential marked as (a), 
TjA = 0.1435 and r\B = 0.2565 for the potential marked as 

(b) , and r/A = 0.35 and t/b = 0.05 for the potential marked as 

(c) . Potentials (b) and (c) are shifted by In |/3W(z) =10 and 
20, respectively. 



+ Ad exp(— a z) cos(fiiz — 6<j) ■ (25) 

In Fig. [7| we show the depletion potentials between 
a single big hard sphere of diameter uq — 3(7^ and a 
semipermeable planar membrane. All potentials were 
evaluated for the constant total packing fraction of the 
smaller spheres, tja +Vb =0.4. The potential marked by 
(a) is for t]a = 0.05 and tjb = 0.35, whereas the potential 
marked by (c) is for r)A = 0.35 and r\B = 0.05. We find 
that for the system (a), rich in B species, the real part 
of the leading order pole a\ — 3.28257er^ 4 1 . This leads 
to the oscillation wavelength 1.91411(7^, which is close 
to <jb- On the other hand, for the system (c), rich in 
A species, a\ = 6.41187CT7 1 and this leads to the oscil- 
lation wavelength 0.979930^1. The asymptotic form [cf. 
Eq. H24fl . dashed lines in Fig.[7j agrees very well with the 
DFT results (solid lines) already for separations of the 
order of several diameters. The depletion potential (b) 
was calculated for tja = 0.1435 and r\s — 0.2565. These 
packing fractions are very close to the crossover point. 
We note that the asymptotic form, Eq. l|25fl. gives an ex- 
cellent approximation already for z > ha a- This is the 
consequence of the fact that the asymptotic form sets 
in at intermediate separations. We conclude that the 
asymptotics of the depletion potential near a semiper- 
meable membrane is the same as for the potential of the 
same mixture near a nonpermeable wall. 

V. CONCLUSIONS 



functions Cy(r) and their Fourier transforms are readily 
obtained from the density functional theory 



i(r) 



-0 



S 2 F ex [p A (r),PB(r)} 
8fH(r)6pj(T') 

^ — } M„w (g) W \ J > 



„ dn a dnp 

a./3 







(23) 



In the above <S> denotes the convolution. Explicit expres- 
sions for Cij(r) can be found in Refi^i. 

The asymptotic decay for pair correlations determines 
the decay of other structural and thermodynamic quan- 
tities describing inhomogeneous fluids. Thus the asymp- 
totics of the depletion potential between a hard-sphere 
colloidal particle and a semipermeable planar membrane, 
for systems away from a crossover point, can be written 
as 



f3W(z) ~ Adexp(— qqz) cos(aiz — Qd) 



(24) 



The amplitude Ad and the phase Qd of the depletion 
potential depend on the external potential whereas the 
characteristic decay length Oq 1 and the oscillation wave- 
length 2n/ai are the same as in the decay of the pair 
correlation functions. However in the proximity of a 
crossover point a two-pole approximation is necessary 

(3W(z) ~ Ad exp(— a z) cos(ai2: — 0<j) 



Within the framework of density functional theory we 
have investigated the effect of the membrane permeabil- 
ity on depletion potentials in colloidal systems. We have 
focused on highly idealized systems and considered hard 
sphere colloids and an infinitely thin ideal membrane. 
The main conclusion is that the colloid-semipermeable 
membrane effective interactions are generally weaker 
than those near a hard nonpermeable wall. For hard- 
sphere-ideal-membrane systems this finding can be as- 
certained by analyzing statistical-mechanical sum rules 
which relate the osmotic pressure to the contact value of 
the density profile of the nonpermeable component. We 
have found that the depletion potential gets stronger with 
an increase of the packing fraction of the nonpermeable 
species and with the osmotic pressure. The structure of 
the depletion potential and its asymptotics is the same as 
for the hard wall systems. In particular, we have found 
the structural crossover behavior that separates two dis- 
tinct asymptotic regimes of the depletion potential. Like- 
wise, the influence of the membrane curvature on the ef- 
fective potential is similar to the hard-wall case but the 
effect is moderated. Although no comparison with com- 
puter simulation is presented we expect our White-Bear 
theory results to be very accurate. 

It follows from the present study that in order to in- 
duce sufficiently strong colloid-membrane attraction the 
membrane should have high enough curvature. This ef- 
fect can be brought about by changing locally the shape 
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of the membrane in the near vicinity of a colloidal particle 
and our results underline the importance of this process. 
Another possibility is to increase the osmotic pressure. 
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